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ABSTRACT 

We present a new theoretical model for the luminosity functions (LFs) of Lya emitting galaxies in the frame- 
work of hierarchical galaxy formation. We extend a semi-analytic model of galaxy formation that reproduces 
a number of observations for local and high-z galaxies, without changing the original model parameters but in- 
troducing a physically-motivated modelUng to describe the escape fraction of Lya photons from host galaxies 
(/esc")- Though a previous study using a hierarchical clustering model simply assumed a constant and universal 
value of ft^", we incorporate two new effects on ft^": extinction by interstellar dust and galaxy-scale out- 
flow induced as a star formation feedback. It is found that the new model nicely reproduces all the observed 
Lya LFs of the Lya emitters (LAEs) at different redshifts in z ^3-6. Especially, the rather surprisingly small 
evolution of the observed LAE Lya LFs compared with the dark halo mass function is naturally reproduced. 
Our model predicts that galaxies with strong outflows and f^^" ^ 1 are dominant in the observed LFs. This 
is also consistent with available observations, while the simple universal fl"^^ model requires /i^c" ^ 1 not to 
overproduce the brightest LAEs. On the other hand, we found that our model significantly overpredicts LAEs 
at z > 6, and absorption of Lya photons by neutral hydrogen in intergalactic medium (IGM) is a reasonable 
interpretation for the discrepancy. This indicates that the IGM neutral fraction jchi rapidly evolves from xhi ^ 1 
at z < 6 to a value of order unity at z ^ 6-7, which is broadly consistent with other observational constraints on 
the reionization history. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — methods: numerical 



1. INTRODUCTION 

In the cold dark matter (CDM) models of structure forma- 
tion, structures are formed and grow hierarchically via grav- 
itational instability; subgalactic clumps are formed in CDM 
halos, and they subsequently merge and collapse to grow into 
more massive galaxies. The ionizing radiation emitted from 
massive stars in primeval galaxies should lead to prominent 
Lya emission via the recombination of hydrogen in their in- 
terstellar medium (ISM). Therefore, as predicted by Partridge 
& Peebles (1967), detecting these Lya emissions via narrow- 
or intermediate-band imaging is one of the most powerful 
technique for seeking high-z young star forming galaxies. In 
the last decade, many Lya emitters (LAEs), which seem to be 
powered by star formation activities, have been found through 
this technique (e.g., Cowie & Hu 1998; Hu et al. 1998; 
Rhoads et al. 2000; Taniguchi et al. 2005; Shimasaku et al. 
2006; lye et al. 2006; Murayama et al. 2007). The Lya lumi- 
nosity functions (LFs) of the LAEs are one of the most funda- 
mental observational quantities and they are becoming more 
firmly confirmed because of the increase of the survey fields 
and available samples Through comparison between the 
observed LFs and theoretical models of the LAEs, we should 
be able to obtain important information for LAEs, and more 
generally, for galaxy formation. 

Nevertheless, theoretical understanding of LAEs is still de- 
veloping, and there are several different approaches for the 
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' In this paper, the term "LAE LP" refers to the Lya line luminosity func- 
tion of LAEs, though some other types of LAE LFs can be defined, e.g., UV 
continuum LF of LAEs. 



theoretical modelling: analytical calculations based on the 
halo mass function and on simple assumptions linking Lya 
emission to halo mass (Haiman & Spaans 1999; Thommes 
& Meisenheimer 2005; Dijkstra et al. 2007b; Mao et al. 
2006; Stark et al. 2007a), cosmological hydrodynamic sim- 
ulations (Barton et al. 2004; Dave et al. 2006; Nagamine et 
al. 2006; Tasitsiomi 2006), and hierarchical clustering mod- 
els of galaxy formation (so-called semi-analytic models; Le 
Delliou et al. 2005, 2006). A key ingredient of the theoretical 
modelling is the escape fraction of the Lya photons from their 
host galaxy, /i^c"- In order to predict it precisely, a detailed 
calculation of radiative transfer of the Lya photons in real- 
istic matter distribution in a galaxy is required (e.g., Hansen 
& Oh 2006). However, such procedure costs huge computa- 
tional time because of the resonant scatterings of Lya photons 
in a galaxy (Neufeld 1990; Chariot & Fall 1993) and a direct 
comparison to the observed LFs is difficult. Therefore, previ- 
ous studies for LAE LFs introduced simple phenomenological 
prescriptions. Le Delliou et al. (2005, 2006) and Dijkstra et 
al. (2007b) assumed a constant and universal escape fraction 
regardless of the physical properties of galaxies. Mao et al. 
(2006) adopted screen-type interstellar dust attenuation, but a 
detailed merger history of dark haloes in hierarchical cluster- 
ing was not taken into account. 

Fortunately, there are some hints for the physical proper- 
ties of LAEs from observations. Lya emitting galaxies in the 
nearby universe generally have low metallicity or small dust 
amount (Chariot & Fall 1993) and show evidence for outflow 
or galactic wind (Lequeux et al. 1995; Kunth et al. 1998, 
2003; Mas-Hesse et al. 2003; Keel 2005). These properties 
are also found in high-z Lyman-break galaxies (LBGs) show- 
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ing strong Lya emission (Pettini et al. 2002; Shapley et al. 
2003; Bower et al. 2004; Wilman et al. 2005; Frye et al. 
2007; Pentericci et al. 2007; Swinbank et al. 2007; Tapken et 
al. 2007). It is theoretically quite reasonable to expect that the 
amount of metal or dust and existence of galaxy-scale outflow 
strongly affect the emergent luminosity of Lya emission. Be- 
cause of the resonant scattering by H I, Lya photons have a 
much longer path to diffuse out from a galaxy than UV con- 
tinuum photons, and hence they are even more vulnerable to 
the absorption by dust. On the other hand, outflows would re- 
sult in a velocity difference between H I gas and Lya photons, 
and hence the escape of a Lya photon would become easier. 
Therefore, these observations and theoretical considerations 
indicate that the universal and constant ftl" in all galaxies is 
oversimplified, and a more realistic modeling of ftlc" is worth 
investigated. 

Another hint for the value of ftlc" comes from direct ob- 
servational estimates. Gawiser et al. (2006) reported observa- 
tional estimates of ftZ"' for LAEs at z = 3. 1 as ftZ"' ^ 0.8 (the 
best fit) and an lower limit of /e';;c" ^ 0.2. However, the semi- 
analytic prediction by Le Delliou et al. (2005, 2006) requires 
a much smaller value of /^c" = 0.02 to fit the observed LAE 
LFs. Therefore a new model of LAE LFs that can reproduce 

the observations with f^^c" ~ 1 is highly desirable. 

Here, we present a new model of LAE Lya LFs, based 
on a semi-analytic model of hierarchical galaxy formation [a 
slightly updated version of Nagashima & Yoshii (2004) used 
in Nagashima et al. (2005)]. This Mitaka model is one of 
the latest semi-analytic models in which the merger history of 
dark matter haloes is taken into account and baryon physics 
such as star formation and feedback is phenomenologically 
treated (e.g., Kauffmann et al. 1993; Cole et al. 1994; Na- 
gashima et al. 1999; Somerville & Primack 1999; Nagashima 
& Yoshii 2004; Baugh et al. 2005; Nagashima et al. 2005). 
This model can reproduce most of the observations for photo- 
metric, kinematic, structural, and chemical properties of local 
galaxies. Moreover, the LFs and angular two-point correlation 
functions of LBGs at z = 4 and 5 predicted by this model cou- 
pled with N-body simulations are found to be in good agree- 
ment with the observation (Kashikawa et al. 2006a). We ex- 
tend this model to treat LAEs, but without changing any orig- 
inal parameters determined by the fit to observations. Rather, 
we try to explain the LAE LFs by minimal extension with 
a least number of new and physically-motivated parameters. 
Specifically, we incorporate two new effects in the calcula- 
tion of /i^c": Lya photon extinction by dust and galactic wind 
driven by the supernova feedback. The new model will be 
carefully compared with all currently available data of high 
redshift LAE LFs. 

The observed number of LAEs would be significantly re- 
duced if the intergalactic medium (IGM) is neutral, because 
the red wing of the Gunn-Peterson trough should attenuate the 
Lya emission (e.g., Miralda-Escude 1998; Haiman & Spaans 
1999). Therefore, the apparent evolution of LAE LF is rec- 
ognized as an invaluable probe of the cosmic reionization, 
and this approach has akeady been applied to the observed 
data (Rhoads & Malhotra 2001; Haiman 2002; Malhotra & 
Rhoads 2004; Stern et al. 2005; Kashikawa et al. 2006; lye 
et al. 2006) to constrain the reionization history, which is 
complementary to the spectra of quasars (Fan et al. 2006) 
or gamma-ray burst (GRB) afterglows (Totani et al. 2006). 
However, a weak point of the LAE-LF method is a degener- 



acy between the intrinsic evolution of the LAE LFs and the 
apparent evolution by the IGM effect. Therefore, the intrinsic 
LF evolution must be known reliably to derive a robust con- 
straint on the reionization. However, an ad hoc assumption of 
no intrinsic evolution was invoked in most of previous stud- 
ies. We will discuss some implications for reionization based 
on our results of LAE LF evolution. 

The rest of this paper is organized as follows. In § |2] we 
describe our extension of the Mitaka model to incorporate 
LAEs, and we compare the model results with the observed 
LAE LFs at various redshifts in3<z<7in§[3] We discuss 
implications for reionization in § HI and then summarize this 
work in § |5] The background cosmology adopted in this pa- 
per is the standard ACDM model: ilu = 0.3, = 0.7, fib = 
0.04, /! = 0.7, and 0-8 = 0.9. 

2. MODEL DESCRIPTION 

The detailed description of the Mitaka semi-analytic model 
for hierarchical galaxy formation is given in Nagashima & 
Yoshii (2004) (see also Nagashima et al. 2005). Here we 
focus on the extension of the original Mitaka model to include 
LAEs. 

2.1. Star Formation Rate in Starburst Galaxies 

In the original Mitaka model, all of the cold gas in a galaxy 
turns into stars and hot gas instantaneously when the galaxy 
is classified as a "starburst" as a result of a major merger, and 
hence the SFR cannot be defined. This means that we cannot 
calculate the Lya photon production rate appropriately that is 
essential to determine the Lya luminosity. However, starburst 
galaxies could have a significant contribution to LAE LFs, 
because they have the highest SFR. Therefore we modify the 
Mitaka model to calculate SFR in starburst populations; we 
adopt an exponential SFR evolution with a time scale of rbmst: 
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Tburst 
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where M^^^^ is the initial mass of available cold gas for star 
formation. We make a reasonable assumption that Tburst is pro- 
portional to the dynamical time of the newly formed spheroid, 
as Tbui-st = /burst Tdyn- We adopt /burst =10, which is consistent 
with SPH simulations (e.g., Kobayashi 2005) and with a re- 
cent observation of submilUmeter galaxies at z ^ 2-3.4 (Tac- 
coni et al. 2006). 

It is expected that not all of the cold gas will be locked 
up into stars, but some fraction of gas will be heated by su- 
pernova feedback, and then ejected from a galaxy as a hot 
galactic wind. This effect has already been taken into account 
in the Mitaka model, to reproduce various scaling laws of lo- 
cal galaxies. We define the efficiency of star formation /* as 
= /*MJ?jj[jj, where is the finally produced stellar mass 
after the starburst episode. The parameter /* is determined 
by the Mitaka model, and it depends on the dark halo circu- 
lar velocity Vc, since the feedback should be more efficient 
for less massive (smaller Vc) galaxies from which gas can be 
more easily removed.^ 

Since the galactic wind will blow up the ISM gas and stop 
the star formation, it is reasonable to identify the time of on- 
set of the galactic wind, fwind, as that when the stellar mass 

2 The quantity /.(V^) scales as <x and oc V° for Vc < K-,. and > y„, 
respectively, where Vc,* ~ 130 km/s (Nagashima & Yoshii 2004). 
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produced by the above SFR evolution becomes f*M'^^^^, i.e.. 



(2) 



and hence fwind = -Tbur,stln(l - /*). This indicates that the 
galactic wind phase will onset later for massive galaxies, and 
this trend is similar to the traditional picture of the galactic 
wind in starbursts (Arimoto & Yoshii 1987). After the on- 
set of the galactic wind, star formation is stopped and the in- 
terstellar gas will escape from a galaxy with a time scale of 
~ fesc = re/Vc, where is the effective radius of the galaxy 
calculated in the Mitaka model. 

2.2. Observed Line Luminosity ofLya 
Let L™" be the Lya line luminosity emitted from a galaxy, 
which is given by 



T emit 
^Lya 



■Lya 
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where ftl,^ and ftZ" are the escape fraction of H I ionizing 
photons and Lya photons from their host galaxy, respectively. 
We adopt flZ'^ = 0, i.e., all ionizing photons emitted from 
massive stars are absorbed in their host galaxy. Although this 
assumption may be rather extreme, it is a reasonable treatment 
for our purpose because the observationally inferred value of 
ftl^ is much less than the unity (Inoue et al. 2006). The mod- 
elling of /^c" is the key issue of this work, and will be treated 
in detail in the next subsection. 

The maximally possible luminosity L™^ is the luminos- 
ity that is achieved if all H I ionizing photons emitted from 
massive stars are absorbed by H I in the ISM (i.e. fl'Z'^ = 0) 
and then reprocessed into Lya photons in the ionization equi- 
librium. We define ryLya as the expected number of Lya 
photons produced by one ionizing photon (or equivalently, 
by one recombination); in the case A recombination it is 
given by ryLya = a'^^/a^ = 0.40 and in the case B, it be- 
comes ryLya = iaB-a%)/aB = 0.68 (Osterbrock 1989, §4.2 
and 11. 8). Here, T = 10'^ K is assumed in the both cases. We 
apply the case B following the previous studies (e.g.. Chariot 
& Fall 1993; Valls-Gabaud 1993; Haiman & Spaans 1999; Le 
DelUou et al. 2005, 2006; Mao et al. 2006; Tasitsiomi 2006; 
Dijksti-a et al. 2007b; Stai'k et al. 2007a), but the luminosity 
would not be changed significantly even if we apply the case 
A. Note that the case A may be appropriate for the situation 
1. 

Then, for a quiescent galaxy in which star formation time 
scale is much longer than the massive star lifetime, L™^ is 
given by 

i^yMt) = ?7LyaeLyaV^(f)er'(Zcold), (4) 

where e^ya = hi'Lya = 10.2 eV and Qh""^^ is the ionizing photon 
emission rate normalized by a unit SFR. On the other hand, 
star formation time scale could be comparable with the mas- 
sive star lifetime in starburst galaxies, and we must take into 
account the evolution of ionizing photon production rate in 
stellar evolution. Then we calculate as 

^i"ya(0=^LyaeLya / V(f')GH Zcold)df', (5) 

Jo 

where 2H""'(^^coid) is the ionizing photon production rate 
from an unit stellar mass whose age is f . Note that the physi- 
cal dimensions of Q'^^'^'' and are different. These quan- 
tities are calculated from the result of Schaerer (2003) with 



of ftic"' 



a correction for the initial mass function used in the Mitaka 
model. Both Qj^""^'' and Qj^""' increase by a factor of 3 with 
decreasing stellar metallicity from the solar abundance to the 
zero metallicity. The stellar metallicity is calculated from that 
of the cold gas, Zcoid, in the Mitaka model. The evolution of 
2y""*'(f ) is characterized by the two phases: almost constant 
phase at 2^"'"=*' ~ 10"^^ photons s"' M©"' for t < 10^ yr and 
the subsequent exponential decay with a typical timescale of 
several million years. 

Finally, when the IGM neutrality is high, Lya luminosity 
could be reduced by IGM absorption, and we include this ef- 
fect as: 

robs _ remitT^IGM 



"^Lya 



"^Lya Lya 



where is the IGM transmission to Lya photons, and it 
is set to be unity except for z > 6, since T^y^ ^ 1 has been 
established at z < 6 by observations (Fan et al. 2006). 

2.3. Lya Escape Fraction Modelling 
We first examine the following models for ftlc'- 

• the simply proportional model: ftZ"' = /o, 

• the dust (screen) model: fs^c"' = /oe"'^. 



the dust (slab) model: /e 



•Lya 



:/o(l-e-')/x. 



where x = A^coidZcoid/ (A^coidZcoid)o. and A^coid is column den- 
sity of cold gas that is computed from the cold gas mass Mcoid 
and the effective radius r^. via A'coid = (■Wcoid/2) /Trr^. The 
simply proportional model assumes an universal and constant 
value of fs^c" for all galaxies, as assumed in Le Delliou et 
al. (2005, 2006). The parameter (A'coidZcoid){) determines the 
strength of the absorption of Lya by dust. The absorption of 
UV continuum by dust has already been included in the Mi- 
taka model, but we treat this parameter independently for the 
Lya absorption by dust, since the effective optical depth for 
the Lya photons could be much larger than that for contin- 
uum photons because of much longer path caused by multiple 
resonant scattering of Lya photons by neutral hydrogen. 

Furthermore, we examine another /i^c" model that incorpo- 
rates the effect of galaxy-scale outflow in addition to that of 
the interstellar dust extinction, which we call the outflow+dust 
model. In this model, starburst galaxies are classified into 
three phases of pre-outflow, outflow, and post-outflow phases. 
These phases are defined by the elapsed time t from the on- 
set of star formation, as t < fwind, fwind < ? < ?wind + fesc, and 
f > fwind + fesc, respectively (see § 2.1 for the definition of 
fwind and fesc)- All the quiescently star-forming galaxies are 
classified into the pre-outflow phase. We model fl 
outflow+dust model as follows: 

• the pre-outflow phase: /i^^c" = fosT'', 

• the outflow phase: ftZ" = f^""^. 



■Lya 



in the 



the post-outflow phase: Lf^^ = 



In the pre-outflow phase, we adopt the same modeling as the 
above dust models. We only apply the screen dust prescrip- 
tion because the difference between the screen and slab dust 
models is small, as we will find in § |3] In the outflow phase, 
according to the theoretical expectation that outflow could 
drastically decrease the effective opacity for Lya photons, we 
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adopt /(*'" = 0.8 as inferred for observed LAEs (Gawiser et 
al. 2006). Note that we do not treat /q""^ as a free parameter, 
keeping the number of free parameters same as the dust mod- 
els. In the post-outflow phase, while the H I ionizing photons 
are produced at moderate rates, neutral gas to absorb these 
photons in host galaxies is absent (i.e., ftl,^ = 1), and hence 
no Lya photons are produced (L™" = 0). 

Therefore, the new free parameters in all the above models 
are /o and (A^coid-Zcoid)o^ while all other parameters are from 
the original Mitaka model without changing their values. 

2.4. Luminosity Function of LAEs 

Now we can calculate the Lya LF of LAEs by calculat- 
ing LLya explained above from the Mitaka model of galaxy 

formation. We determine the free parameters of ftsc"' mod- 
els by fitting the model Lya LF with T^^^ = 1 to the data 
at z = 5.7 reported by Shimasaku et al. (2006), which is cor- 
rected for both detection completeness and contamination, be- 
ing the most reliable LAE LF to date. Specifically, we per- 
form a test in the whole observational range of the Lya 
luminosity (log [Li,ya/(h~^ ergs s"')] =42.1-43.1). Because 
the LF error bars are the smallest at the faint end, the free pa- 
rameters of ftsc' are determined mainly by matching to the 
faint-end of the Lya LF of Shimasaku et al. We fix these val- 
ues of the free parameters at z = 5.7 and then apply the model 
to the data at other redshifts with keeping the parameter val- 
ues, because they should reflect the physical properties of Lya 
photons that are not expected to evolve with redshift. 

It should be noted here that we do not set any criterion on 
the equivalent width (EW) of Lya lines for model galaxies 
to be selected as LAEs. Though the observations of LAEs 
usually set threshold values for EW in the selection, the ef- 
fect of different EW thresholds on the predicted Lya LF is 
small. We show the model calculation of Lya EW distribu- 
tion of LAEs at z = 5.7 in Figure [T] which is intrinsic EW 
and the IGM absorption effect on the continuum is not taken 
into account. Here, we calculated the EW by L^^a ^^'^ '^^e 
continuum luminosity of model galaxies. Almost all LAEs 
have larger EWs than a typical threshold in observation, i.e., 
EW^'^\Lya) = 20 A, and hence our results do not change even 
if we set a threshold Lya EW used in observations. This result 
is consistent with that of Le Delliou et al. (2006). 

3. RESULTS 

3.1. Luminosity Functions at 6 

We first compare the simplest model (the simply propor- 
tional model) to the cumulative Lya LF data at z = 5.7 in Fig- 
ure |2] The best-fit is /o = 0.60, and this model overproduces 
the bright-end of the LAE LF. We see from the top panel of 
Figure|2]that the dominant component at the bright-end of the 
Lya LF is the starburst galaxies, while the quiescently star 
forming galaxies dominate at the faint end. If we try to adjust 
the model to the bright end, we must assume a much lower 

value of ftlc" ^ 1 than the observationally inferred value by 
Gawiser et al. (2006). This could be one of the possible rea- 
sons that Le Delliou et al. (2005, 2006) obtained fts" < 1 
with the same treatment for fsZ"' , though their model is com- 
pletely independent of our own and a clear comparison is not 
easy. 

The middle and bottom panels show that the starburst galax- 
ies in the pre-outflow phase having high metal column densi- 



ties (A^coid-Zcoid > 10^' Zq cm"^) dominate in the largest lumi- 
nosity range. This is reasonable since the most efficient Lya 
photon production is expected during young stage of the star- 
burst, and the galaxies in such phase are expected to be gas 
and metal rich. This result indicates that the effect of interstel- 
lar dust extinction in ftZ"' models would reduce the number 
of LAEs at the bright end, in a direction to a better agreement 
with observations. 

We then show the dust models and the outflow+dust model 
in Figure [3] The fitted values of the model parameters are 
tabulated in Table [T] As expected, the number of the lumi- 
nous LAEs drastically decreases in the dust models. However, 
it can be seen that this effect is too strong, and the sharp drop 
of LAE number at the bright end does not match well to the 
observed Lya LFs. In contrast, the outflow+dust model well 
reproduces the observed Lya LFs in the whole range of Lya 
luminosity shown in Figure [3] It should be noted again that 
the number of free model parameters is the same for the dust 
and outflow+dust models. LAEs at the bright-end are dom- 
inated by galaxies in the outflow phase, which is consistent 
with the observational indications of outflows. Furthermore, 
/eiic" 0.8 has been assumed for the outflowing LAEs, which 
is consistent with a recent observational estimate for LAEs at 
z ~ 3.1 by Gawiser et al. (2006). 

3.2. Luminosity Functions at z^6 and IGM transmission 

Turning now to the Lya LFs at redshift beyond 6 when 
the reionization is believed not to have ended yet. In con- 
trast to the agreement at the lower redshifts, all /i^c" models 
assuming Ij^™ = 1 with the same parameters determined at 
z = 5.7 overpredict the number of the LAEs at z = 6.56 com- 
pared with the observational data reported by Kashikawa et al. 
(2006b), as presented in the top panel of Figure |4] This dis- 
crepancy can be resolved if we adopt a simple prescription of 
luminosity-independent IGM transmission that is less than the 
unity (Tl^^ - 0.6-0.8); the predicted Lya LFs with T^^^ = 
0.8 and 0.6 are shown in the middle and bottom panels, re- 
spectively, in Figure m We also show the predicted Lya LFs 
at an even higher redshift of z = 6.96 in Figure |5] Though 
the statistics is obviously insufficient, the most distant galaxy 
found so far (lye et al. 2006) indicates 7]^™ < 0.7. We will 
discuss the implications for the reionization history from these 
results in § |4] 

3.3. Interpretations for Intrinsic Lya LF Evolutionary 
Features 

We show the intrinsic (i.e., T^y^ =1) evolution of the cumu- 
lative Lya LF in the outflow+dust model from redshift z = 3.6 
to 8.8 in Figure|6] It is also interesting to compare this evolu- 
tion with that of the dark halo mass function, since a num- 
ber of previous studies predicted LAE LF evolution based 
on the halo mass function. This is shown in Fig. |7] using 
the halo mass function given by Yahagi et al. (2004, here- 
after YNY) that is used in the Mitaka model. Here, we adopt 
MoM/LLya = 160 Mq/Lq, whcrc Lq is the bolometric solar 
luminosity. 

A clear trend of the LAE LF evolution is that the charac- 
teristic luminosity where the LFs have a break, Lbr (similar to 
in the Schechter function), shows very little evolution in a 
wide range of redshift. This seems the primary reason that the 
LF is not very sensitive to the redshift at Llyq ^ L^r- On the 
other hand, the degree of the evolution at LLya ^ ^br is similar 
to that of the dark halo mass function. Therefore, a mecha- 
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nism that keeps the brightest Lya luminosity almost constant 
against redshift would explain these trends. 

Such a mechanism is likely to be a combination of the dust 
extinction and outflow effect. The comparison of the dust 
and dust+outflow models with the simply proportional model 
(Figs. |2]and[3]i indicates that galaxies more massive than a 
critical mass scale cannot become bright LAEs because of 
heavy extinction by a large amount of gas and metals, and/or 
because of the large gravitational potential that prevents an 
efficient outflow. Therefore, even if the number of galaxies 
more massive than the critical mass scale should increase sig- 
nificantly with time by hierarchical structure formation, the 
mass of the brightest LAEs does not significantly evolve. This 
interpretation can be tested by estimating the mass of LAEs 
in future observations. 

4. IMPLICATIONS FOR COSMIC REIONIZATION 

The requirement of Tj^™ < 1 to reproduce the observed 
Lya LFs at z > 6 suggests that the IGM opacity for the Lya 
emission rapidly increases beyond z ^ 6. From the value of 
^Ly'a' estimated in our analyses, we can estimate the IGM neu- 
tral fraction xm = «Hi/nH- Unfortunately, this procedure is 
highly complicated because there are many physical processes 
that must be considered in calculation of the attenuation factor 
of Lya line luminosity (Santos 2004; Dijkstra et al. 2007a). 
Here, we apply the dynamic model with a reasonable velocity 
shift of Lya line by 360 km/s redward of the systemic velocity 
(the dashed curve of Figure 25 in Santos 2004). 

The reason for the choice of this particular model is that it 
predicts no attenuation when xhi = 0. Note that some other 
models of Santos (2004) predict a significant attenuation even 
in the case of xui = 0, due to the neutral gas associated with 
the host haloes of LAEs. Choosing this model then means 
that we ascribe the sudden strong evolution of the Lya LF at 
z > 6 only to the absorption by pure IGM. We consider that 
this is a reasonable assumption, since observations indicate 
that the escape fraction of Lya photons is about unity at least 
for LAEs at z - 3 (Gawiser et al. 2006). If LAEs at z - 7 
are a similar population to the low-z LAEs, we do not expect 
significant absorption by neutral gas physically associated to 
LAEs. It should also be noted that the brightest LAEs in our 
model LFs are in the outflow phase, and we may not expect a 
large amount of neutral gas around them. 

For the range of = 0.6-0.8 at z = 6.56 estimated from 
the outflow+dust model, we findxei ~ 0.25-0.35. Atz = 6.9, 
our constraint of 7"^^^^ < 0.7 translates into xhi > 0.3. On the 
other hand, it should also be kept in mind that if z 7 LAEs 
are surrounded by a significant amount of nearby neutral gas 
that is not present for LAEs at z ~ 3, the estimate of xhi as 
an average of IGM in the universe could become lower than 
those derived here. 

Although the translation from T^y^ into xui is model de- 
pendent, it should be noted that our model overpredicts the 
observed LAE LFs rather suddenly beyond z 6, while it re- 
produces well the observed LAE LFs at z < 6 without invok- 
ing the absorption in IGM. We suggest that this is an evidence 
for significant absorption by neutral IGM at z > 6, and it re- 
quires Xhi of order unity. This is broadly consistent with the 
recent constraints by other methods (Fan et al. 2006). Here we 
emphasize that there has been no strong "positive" evidence 
for a considerable amount of IGM neutral hydrogen (xhi of 
order unity) at z > 6; quasar Gunn-Peterson tests give only a 
weak lower limit of xhi > 10"^ (Fan et al. 2006) and the GRB 



050904 gives only an upper limit of xhi < 0.6 (Totani et al. 
2006). 

For the future studies of reionization through LAE LFs, we 
plot the LF evolution of the outflow+dust model from z = 3.6 
to z = 8.8 in Figure|6]with r/°^ = L The redshifts of 7.7 and 
8.8 correspond to the next windows that are relatively free of 
bright OH lines at wavelengths of 1.06 /im and 1.19 fj,m, re- 
spectively. There are already several pioneering works, such 
as Willis & Courbin (2005), Cuby et al. (2007) and Stark et 
al. (2007b). Moreover, other wide-field imagers with narrow- 
band filters to search the redshifted Lya emission at these 
wavelengths will be available on 4 m or 8-10 m class tele- 
scopes in the near future (see, e.g. Cuby et al. 2007, for more 
detail). Our predicted LAE LF evolution would be useful for 
the planning of the future LAE surveys, and it can also be 
used as a guide to derive Tj^™ from the observed LAE LFs^. 

5. CONCLUSIONS 

We constructed a new theoretical model for Lyman alpha 
emitters at high-z based on a hierarchical clustering model 
of galaxy formation (the Mitaka model), taking into account 
physical effects of dust absorption and galaxy-scale outflow 

in calculation of Lya photon escape fraction iftic") from host 
galaxies. We introduced just two new parameters based on 
some physical considerations to include these two effects, 
while we kept unchanged the original model parameters of 
the Mitaka model that have been tuned to reproduce a number 
of observations of local galaxies. 

The observed Lya LFs of LAEs at several redshifts in z < 6 
are reproduced well by our model assuming completely ion- 
ized IGM (the IGM transmission = 1). Especially, it has 
been known that the LAE LFs show very little evolution in 
z ^ 3-6 compared with the dark halo mass function, and our 
model can naturally explain this by the above effects newly 
included in this work. The escape fraction of Lya photons 

from LAEs in our model is ftZ°' ^ 0.8, which is consistent 
with a recent observational estimate. Moreover, our model 
predicts that the bright-end of the Lya LF is dominated by the 
starburst galaxies in the phase of galaxy-scale outflow driven 
by the supernova feedback. This is also consistent with the 
observational results for local and high-z galaxies with strong 
Lya emission. On the other hand, the simple model with a 
constant and universal ftZ"' requires f^^^" ^ 1 in order not 
to overproduce the brightest LAEs. This is consistent with a 
previous theoretical study that invoked the same assumption, 
but the required ftZ"' value is considerably lower than the ob- 
served values. 

In contrast to the success at z < 6, our model overpredicts 
the LAE LF beyond z > 6, and we interpret this as a result of 
Lya attenuation by neutral hydrogen in IGM. We found that 
a simple prescription of luminosity-independent Tj^™ is suf- 
ficient to explain the observed LAE LFs; the LF at z = 6.56 
is reproduced with = 0.6-0.8, and that at z = 6.9 indi- 
cates T^y^ ^ 0.7. Though it is not straightforward to derive 
the IGM neutral fraction xhi from these results, the discrep- 
ancy in the case of = 1 appears rather suddenly beyond 
z > 6 compared with the results in z < 6. Therefore we sug- 
gest that this gives an evidence for a rapid evolution of xhi 
from Xhi ^1 at z < 6 to a value of order unity at z > 6. This 

' The numerical data of the LF evolution is available on request to the 
authors. 
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positive evidence for a significant abundance of neutral hy- 
drogen in IGM is complementary to the constraints obtained 
by other methods based on quasar and GRB spectra. 

Our theoretical model for the LAE LF evolution would be 
helpful for planning a LAE survey and for interpretation of the 
observed Lya LF at even higher redshifts in the future stud- 
ies. Further investigation of various properties (e.g., spectral 
energy distribution, mass, or spatial clustering) of LAEs in 
our model and comparison with observations will be done in 
our future studies. 
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TABLE 1 

The best-fit model parameters of Lyo photon escape fraction. 



flZ" model fo (A'coidZcoid)o [lO^° Zq cm ^] 

simply proportional 0.60;!^ S^ 

dust (screen) ^■^''^^ 1 ^■'^-2 

dust (slab) 

outflow+dust 0.68!^-S^ 4:2+\i 



Note. — The errors are statistical Icr as a result of the tit to the Lya LF of the LAEs at z = 5.7 reported by Shimasaku et al. (2006). The screen dust extinction is adopted in the 
pre-outilow phase of the outflow+dust model. 
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Fig. 1 . — The predicted rest-frame EW of the Lyo emission line for the LAEs at z = 5.7 in the outflow+dust model, decomposed into three intervals of Lyo 
luminosity. The EWs are intrinsic without IGM absolution for Lya or UV continuum flux. The thin and thick histograms are the EW distributions for quiescent 
and starburst galaxies, respectively. The vertical dotted line represents a typical threshold for galaxies to be selected as LAEs in observation. 
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log[LL /(h-^erg/s)] 
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log[LLy„/(IT2 erg/s)] 

Fig. 2. — The cumulative Lyo LF of the simply proportional model at z = 5.7 with /o = 0.60. Top: The contiibutions from quiescent and starburst galaxies are 
shown separately, as indicated in the figure. The symbols with error bars are observational data of the LAEs at z = 5.7. The references for the data points are as 
follows: Shimasaku et al. (2006), Murayama et al. (2007), Ajiki et al. (2003), Hu et al. (2004), and Rhoads et al. (2003). Note that these observational data are 
not corrected for detection completeness and contamination, except for Shimasaku et al. (2006). Middle: starburst (thick Hues) and quiescent (thin lines) galaxies 
are further decomposed into two groups corresponding to high (solid lines) and low (dashed lines) values of the metal column density, WcoidZcoldi separated by 
^cold^cold = 10^' [Zq cm"^]. The dotted curve is the total Lya LF without any decomposition. Bottom: The Lya LFs of starburst galaxies are shown by thick 
curves, which are decomposed into the pre-outflow, outflow, and post-outflow phases. The Lya LF of the quiescent galaxies is shown by the thin dashed line, and 
the thin dotted line is the total Lya LF. 
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log[LLy„/(h-^ erg/s)] 

Fig. 3. — The evolution of the cumulative Ly« LFs with redshift at z < 6. The curves show the model results, while the symbols with en'or bars are the 
observational data. The four different models of f^^" , i.e., the simply proportional, dust (screen), dust (slab), and outflow+dust models are shown (see the top 
panel for the line mai'kings). The thin solid curve is the contribution from the starbursts in the outflow phase in the outflow+dust model. The references for the 
data points that were not given in Figure|2]are: Yamada et al. (2005), Gronwall et al. (2007), Gawiser et al. (2006), Kudritzki et al. (2000), Cowie & Hu (1998), 
Maier et al. (2003), Ouchi et al. (2003), and Santos et al. (2004). 
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Fig. 4. — The cumulative Lya LFs at z = 6.56. The hne-markings are the same as Fig. [5] but with different values of the IGM transmission 7^^^^ = 1.0, 0.8, 
and 0.6 for the top, middle, and bottom panels, respectively. The open squares are the observed data from the photometric sample of Kashikawa et al. (2006b) 
corrected for detection completeness (an upper limit), while the circles from the pure spectroscopic sample (a lower limit). 
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log[LLy^/(h-^erg/s)] 

Fig. 6. — The evolution of the cumulative Lya LF in the outflow+dust model from redshift z = 3.6 to z = 8.8, assuming transparent IGM to Lya photons 
(T^'^^ = 1). The contributions from starburst and quiescent galaxies are presented separately in the bottom panel. 
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Fig. 7. — A comparison of the evolution of the cumulative Lya LF in the outflow+dust model (thin curves) with that of the dark halo mass function (thick 
curves) from redshift z = 3.6 to z = 6.9. The conversion factor of the halo mass into the hya line luminosity, Mou/LLya = 160 Mq/Lq, is chosen to fit the 
cumulative hya LF at z = 6.9. 



